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False  Sidelobe  and  Multiple  Target 
Resolution  Probiem 
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Adaptive  beamforming/filtering  and  STAP  methods  can  result  in 
high  sidelobe  levels,  especially  with  limited  sample  support  for 
estimating  interference  covariance  matrix 


This  can  cause  excessive  “false”  sidelobe  detections  arising  from 
targets  or  undernulled  interferences 


Sidelobe  rejection  capabilities  of  adaptive  matched  filter  (AMF), 

generalized  likelihood  ratio  test  (GLRT)  and  AMF/ACE  have  been 
previously  analyzed 


A.  Jaffer  4/9/01 


Multiple  Target  and  Sidelobe  Detection 
using  Weighted  Least-Squares  Fit  to 
AMF  and  GLRT  Data 


Raytheon 


For  simplicity  of  exposition,  we  will  consider  the  spatial  domain 
only  although  the  method  is  directly  applicable  to  the  angle- 
Doppler  domain 

The  proposed  method  uses  the  output  detection  test  statistic 
computed  over  the  entire  angular  extent  of  interest 

In  any  given  range  cell  the  question  is  whether  the  totality  of 
output  test  statisitic  values  computed  over  the  beam  directions  of 
interest  that  exceed  a  preset  threshold  represent  one  target,  two 
targets  or  up  to  a  maximum  M  targets,  i.e.,  the  question  is  one  of 
model  order  determination 
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J-amf  and  J-gIrt  Functions  -  Normalized 
Adaptive  Array  Outputs 
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J-amf(6 )  represents  the  normalized  adaptive  array  power  output 
as  a  function  of  the  angie  0 

H’"(0)RM'(e) 

J-glrt(6)  is  J-amf(6 )  further  normalized  by  data  power  term 


d’‘(Q)R-'x 


GLRT 


(0)= 


AMF 


(e) 


l+x"R~^x/K 


-  1  ^  H 

R  =  —  X}^  ’  s  are  secondary  set  of  signal  -  free  data 

^  k=i 


K  -  -  Number  of  snapshots 
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Normalized  Adaptive  Array  Outputs 

—  continued 
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Diagonal  Loading  and  shaded  (tapered)  steering  vector: 


R[,l  —R  +  oI 


=w  •  represents  the  Schur  product ,  w  is  a  shading  or  taper  function 


w(Q) = Rold^i,  is  also  used  in  the  J-amf(6 )  and  J-glrt(9 )  functions 
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J-amf  Function  for  One  and  Two  Targets 
---No Taper,  K=100 
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Threshold  for 
Pfa=10E(-6) 


Threshold  for 
Pfa=10E(-6) 
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J-amf  function  -  unshaded  steering  vector 
and  50  dB  Chebyshev  taper 


J^mf  (®)  source  at  0  degree,  jammer  at  -30  degree,  using  Rhat  =  sample-average  R,  K  =  1 00. 


Threshold  for 
Pfa=  10E(-6) 


•  False  sidelobe  detections  can  occur  in  adaptive  arrays  even  with 
50  dB  Chebyshev  taper 
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Weighted  Least-Squares  and  Model  Order 

Determination  Method 


Raytheon 


*  Determine  L  peaks  ( or  values  adjacent  to  peaks)  of  Jamf  (o) 
at  61,62, . ,62.  which  exceed  preset  threshold  for  given 


•  Test  data  vector  x  is  modeled  as 


x  =  D^a  +n 


where  ),</(6,2), . ,</(6j„  )J  an  N  by  M  matrix 

a  and  "  are  the  complex  signal  amplitude  vector  and  interference 
plus  noise  vectors 

*  The  application  of  the  weight  vector  >^(6)  to  the  data  yields 

y( 6)  =  X  =  wCS)*  D^a+v 


A.  Jaffer  4/9/01 


WLS  Method . continued 
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where  v(0)=m'(0)®« 


is  evaluated  at  L  distinct  points 

The  transformed  signal  model  is  fit  to  the  data  in  a  weighted 
least-squares  sense  and  the  residual  is  evaluated 

fjil  rM»(0,  fvfOi)! 

I  .  M  .  I  I  .  I 

I  LI  l.l  I 

I  ■  I  I  •  I  I  ■  I 

L'z.J  l_v(0i)J 

The  residual  is  computed  as 

J res  ||^  )|| 


or  Y=Ha+V 
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WLS  Method 


continued2 
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)  is  an  orthogonal  projection  operator 
For  details,  see  paper 

^  J  res  (®sl  V  •  ®sM  )  is  a  function  of  and  needs  to  be  minimized 

over  those  angles 


For  computer  simulation  purposes,  we  limit  ourselves  to  M=2 

As  a  first  approximation,  we  take  o.i  and  6,2  to  be  the  two  highest 
peaks  —  computationally  efficient 

Better  approximation  --  fix  one  angle  at  global  peak  and  search 
over  other  to  minimize  —  more  computation 
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Akaike  Information  Criterion  (AlC)  for 
Model  Order  Determination 
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Compute  AlC(m)  for  model  order  m=1 ,2 . M  and  choose  the 

minimum 


AiC(m)=  Jres  +  Number  of  free  parameters  in  model 
=  Jres*  3m 


Alternate  Minimum  Description  Length  (MDL)  criterion  is  not 
applicable  here  because  the  second  term  ( the  “penalty”  term) 
becomes  zero  for  single  target  signal  snapshot 
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Performance  Evaluation  by  Monte-Carlo 
Computer  Simulation 
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10  element  line  array  with  half-wavelength  inter -element  spacing 

Single  source  placed  at  0  degrees  (broadside) 

Two  sources  at  0  and  45  degrees 

Noise  jammer  placed  at  -30  degrees,  JNR=40dB 

Thresholds  for  AMF  and  GLRT  methods  to  yield  a  specified 
computed  in  accordance  with  paper  by  F.C.  Robey  et.al., 

IEEE  Trans,  on  AES,  March  1992 

The  probability  of  detection  was  based  on  5000  trials  for  each 
point  on  the  curve 

A  target  detection  was  considered  vaiid  if  it  feil  within  pius  or 
minus  3  dB  beamwidth  of  the  true  target  angle 
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PROBABILITY  OF  DETECTION 


Probability  of  Mainlobe  Detection  (AMF) 


AMF:  SINGLE  SOURCE  DETECTION,  N  =  10,  =  10'® 

FA 
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PROBABILITY  OF  FALSE  SIDELOBE  DETECTIONS 


Probability  of  Sidelobe  Detections  (AMF) 


AMF:  FALSE  SIDELOBE  DETECTIONS,  N  =  10,  Pp^  =  10‘® 


A.  Jaffer  4/9/01 


PROBABILITY  OF  DETECTION 


Probability  of  Mainlobe  Detection 
—  AMF  Taper 
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AMF:  SINGLE  SOURCE  DETECTION,  N  =  10,  Pp^=  10'® 
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PROBABILITY  OF  FALSE  SIDELOBE  DETECTIONS 


Probability  of  Sidelobe  Detection 
—  AMF  Taper 


AMF:  FALSE  SIDELOBE  DETECTIONS,  N  =  1 0,  P  =10'® 

FA 
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Probability  of  Mainlobe  Detection 
AMF  Taper  plus  Diagonal  Loading 


AMF:  SINGLE  SOURCE  DETECTION  (DIAGONAL  LOADING  +  TAPER),  N  =  10,  K  =  20,  =  10'® 
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Probabilty  of  Sidelobe  Detection 
AMF  Taper  plus  Diagonal  Loading 


AMF:  FALSE  SIDELOBE  DETECTIONS  (DIAGONAL  LOADING  +  TAPER),  N  =  10,  K  =  20,  Pp^=  10'^ 
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PROBABILITY  OF  DETECTIONS 


Probability  of  Detecting  Two  Sources  -  AMF 


AMF:  TWO  SOURCES  DETECTIONS,  N  =  10,  =  10'® 
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PROBABILITY  OF  DETECTIONS 


Probability  of  Detecting  Two  sources  -  GLRT 


GLRT:  TWO  SOURCES  DETECTIONS,  N  =  1 0,  P  =10'® 

FA 
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PROBABILITY  OF  DETECTION 


Probability  of  Detection  -  -  GLRT 


GLRT:  SINGLE  SOURCE  DETECTION,  N  =  1 0,  P  =10'® 

FA 
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PROBABILITY  OF  FALSE  SIDELOBE  DETECTIONS 


Probability  of  Sidelobe  Detection 


GLRT 


GLRT:  FALSE  SIDELOBE  DETECTIONS,  N  =  10,  P  =10'® 

FA 


A.  Jaffer  4/9/01 


PROBABILITY  OF  DETECTIONS 


Probability  of  Detecting  Two  Sources 
-  -  -  Various  MT-AMF  methods 
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AMF:  TWO  SOURCES  DETECTIONS,  N  =  32,  K  =  64,  Pp^=10'® 
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ARRAY  SNR  OF  THE  SIDELOBE  SOURCE  (dB) 
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Proposed  new  method,  using  weighted  least-squares  fit  to  AMF 
or  GLRT  data  combined  with  model  order  determination  by  Akaike 
Information  Criterion,  can  significantly  reduce  false  sidelobe 
detections 

This  is  true  even  when  amplitude  taper  and  diagonal  loading  is 
used 

The  approximate  method  for  angle  estimation  is  computationally 
efficient  and  yields  good  detection  and  sidelobe  rejection  results 
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